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Abstract. We investigate the tricritical Ising model in complex magnetic field in 
order to characterize the analytic structure of its free energy. By supplementing 
analytic methods with the truncation of conformal space technique we obtain 
nonperturbative data even if the field theories we consider are not integrablc. 
The existence of edge singularities analogous to the Yang-Lee points in the Ising 
field theory is confirmed. A surprising result, due to the conformal dimensions 
of the operators involved, is the appearance of two branching points which seems 
appealing to identify with a pair of complex conjugate spinodal singularities. 
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1. Introduction 

More than fifty years ago, Yang and Lee [f ,2] showed the importance of understanding 
the analytic structure of the free energy in the theory of phase transitions. They 
proved that the thermodynamical equation of state is completely determined by the 
distribution of roots of the partition function (or, equivalently, singularities of the free 
energy). Their papers inspired several lines of further research, which we will not 
attempt to review here. We just note that, as far as the Ising model is concerned, the 
task of unveiling the analytic properties of the free energy has been recently completed 
by Fonseca and Alexander Zamolodchikov [3] . They have been able to substantiate the 
identification of the Yang-Lee edge singularity with the spinodal point, thus unifying 
the high- and low-temperature descriptions. 

One of the main tools of Fonseca and Zamolodchikov's analysis is the truncated 
free-fermion space approach, a technique that gives numerical access to the lowest 
energy levels of the Ising field theory defined on a cylinder. A similar analysis can, in 
principle, be performed to study the next-to-simplest unitary minimal model, namely 
the tricritical Ising model (in the following often referred to by the acronym TIM). Due 
to the universality principle, the insight gathered by working on this conformal field 
theory and its perturbations applies to all tricritical phenomena in two dimensional 
space, characterized by the same Z2 symmetry of the order parameter. 
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Figure 1. Ising model: branch cuts of the free energy in the complex _ff-plane. 
The magnetic field is measured in units of 1 27rc?2 1 15 / s . In the left panel, the 
high-temperature regime is considered and the Yang-Lee edge singularities are 
depicted. The right panel shows the Langer branch cut in the low-temperature 
regime. 

In practice, however, the extension of Fonseca-Zamolodchikov's method is far 
from trivial because the tricritical Ising model possesses four primary fields that are 
relevant in the renormalization group sense, instead of two like in the Ising case. On 
technical ground, moreover, it is not possible to exploit the equivalence between the 
magnetically perturbed Ising theory and a theory of free fermions, so that one is forced 
to employ a less powerful approximation, the truncation of the conformal space of the 
minimal model A^(5,4). Nonetheless, some picture can be outlined, and the resulting 
sketch constitutes the subject of this article. 

In order to make the paper reasonably self-contained, the main points of Fonseca 
and Zamolodchikov's work, as well as earlier progress on the analytic structure of 
the Ising model in complex magnetic field, are collected in the next subsection. 
Section [2] briefly presents the truncation of conformal space (TCS for the sake of 
brevity) technique that we adopt, while the most relevant known facts about the 
tricritical Ising model are summarized in section [3] The results of our investigations 
are thoroughly reported in section U and then synthesized and commented in the final 
section. 

1.1. The Ising model in complex magnetic field 

As written in the introduction, the study of the Ising model in a complex magnetic 
field H was introduced by Yang and Lee in 1952. In their first paper [1] they studied 
the grand partition function of a lattice gas with complex fugacity y, proving that the 
equation of state can be deduced by the distribution of the singularities of the free 
energy F. The second paper [2] established the equivalence between the Ising model 
and a lattice gas with a suitable potential. In particular, the fugacity was found to be 
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Figure 2. Magnetization M versus the magnetic field H for /3 > /3 C . The red 
solid line indicates the stable branches, the green dashed line are the metastable 
branches, while the blue dotted line represents the unstable branch. The value 
H = H Bp where the metastable branch turns into unstable is called spinodal 
singularity. 



proportional to the exponential of the magnetic field: 

yoce" 2 ^. (1) 

Then a theorem was proved that, for a broad class of potentials (including the one 
corresponding to the Ising model), all the zeros of the grand partition function lie on 
the unit circle of the complex y-plane (equivalently, all the singularities of the free 
energy are located on the imaginary axis of the complex i7-plane). 

For f3 < (3 C (the high-temperature regime), the logarithmic singularities of the 
free energy accumulate in the thermodynamic limit towards the points H = ±iHo(f3) 
with real Hq > 0. To visualize the situation we can join these accumulation points by 



drawing a cut along the imaginary axis (see figure 1(a) ). The cut passes through the 
point at infinity: the free energy is analytic for any real value of H because we know 
there is no phase transition. 

The gap between the branch points reduces while lowering the temperature, 
until, at the critical temperature, it happens that Hq((3 c ) = 0. Now it is no more 
possible to analytically continue F(f3, H) from positive to negative H : this is the 
signal of the phase transition. Fisher [4] named Yang-Lee edge singularities the 
points ±iiJo(/3), and proved that the accumulation of zeros of the partition function 
in the thermodynamic limit is a conventional critical phenomenon, with scaling laws, 
universality and so on. In high dimension D, it corresponds to the infrared behaviour 
of the field theory with an action 

^(d^) 2 + i(h - h )(j> + i 7 : 



Ayl = I d D x 



(2) 



Note that such a theory is nonunitary because of the imaginary couplings. The theory 
of Yang-Lee edge singularity was enriched by Cardy [5] who showed how, at the critical 
point, it is related to the simplest nonunitary minimal model characterized by central 
charge c = —22/5. Away from the critical point, the Yang-Lee edge singularity is an 
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integrable theory, whose exact 5-matrix was proposed by Cardy and Mussardo [6] and 
its Form Factors and correlation functions were computed by Al. Zamolodochikov [7] . 

On the other hand, for (3 > (3 C (low-temperature regime), if one tries to 
continuously change the magnetic field from H > to H < 0, at H — one enters a 
metastable phase (that is, the system is in a local free energy minimum). Eventually 
the spinodal point is reached, where the metastable phase becomes unstable (see 
figure [2]). This point is a singularity, in the sense that (at least in a mean field 
approximation) the isothermal susceptibility diverges as one approaches it. In fact, 
the above picture must be corrected by taking into account thermic fluctuations [8]. 
They make the system start decaying through nucleation before reaching the spinodal 
singularity. A cut (named Langer's branch cut) is therefore drawn along the negative 



H-axis (see figure 1(b)), starting from H — where is a weak singularity (see [9] for 
a mean field theory of spinodal points and first order phase transitions). 

What happens to the spinodal point when the Langer's branch cut is opened? 
Fonseca and Zamolodchikov answer that it is pushed under the cut in order to reappear 
in the high-temperature regime under the name of Yang-Lee edge singularity. This 
identification completes!!! the task of describing the analytical properties of the free 
energy in the complex magnetic field plane. 

In order to support their claim, Fonseca and Zamolodchikov use the Ising field 
theory, formally defined by the action 



A 



FF 



Aift = A( c=1 / 2 ) +9i j a(x)d x + g 2 j e(x)d x. (3) 

In the above formula -4( c=1 / 2 ) represents the action of the simplest unitary minimal 
model A4(4, 3), a(x) and e{x) are the primary fields of the theory, that we can interpret 
as spin and thermal operators, while g± and g 2 are couplings related to departures 
from the critical point in the magnetic or in the temperature direction, respectively. 
Actually, it is possible to put together the first and the last terms on the RHS of ([3]) 
and substitute them with the free fermions action 

— J [ip&ip + 4>d4> + im^V] d 2 a; > (4) 

where m = 2Trg 2 , and ip,ip are the two chiral components of a Majorana field. The 
operator a admits as well a representation in terms of the Majorana fermion and it 
is easy to compute its matrix elements - at a finite volume - on the multi-particle 
fermionic states [3, 12]. Whichever formulation one employs, when the Ising field 
theory is put on a cylinder, its radius R acts as an infrared regulator, so that the 
spectrum is infrared finite. By means of a suitable truncation of the Hilbcrt space the 
problem of measuring the energy levels is reduced to the diagonalization of a finite 
dimensional Hamiltonian. One has access to the free energy by studying the infrared 
(that is, R — > oo) behaviour of the ground state. 



2. Truncation of the conformal space 

The truncation of conformal space approach has been invented by Yurov and Alexei 
Zamolodchikov [13]. It is an approximate method that gives numerical access to the 
spectrum of two-dimensional off-critical theories arising from deformation of minimal 

§ It seems that the first to speculate about a possible connection between spinodal points and the 
Yang-Lee edge singularity was Klein in [10] . 

|| In fact, the comparison with lattice results (see [11]) could still hide some tricky point. 




Figure 3. Spectrum of the theory .AJ : the first 12 levels. 
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Figure 4. Spectrum of the theory .Aj : the black lines are the exact masses, the 
colored ones are energy levels differences AEi = Ej — Eq as computed by TCS. 



models. Unlike other methods, however, the TCS works equally well for integrable 
and non-integrable deformations. 

The Hamiltonian acting on the Hilbert space is of course infinite dimensional; we 
need to truncate the conformal space if we want to be able to compute anything. In 
this paper the truncation is performed by discarding the states of level 6 or more in 
each Verma module. For the tricritical Ising model, this amounts to keep 228 states. 
The suitable truncated basis of the Hilbert space and the matrix elements of Hq and 
Vi have been computed by means of a program for Mathematica described in [14]. The 
numerical diagonalization of the Hamiltonians was performed by the NAG routines of 
Maple. 

This is a good place to say something about the precision of this approach. By 
following a standard procedure, introduced in [15] and recently theoretically justified 
[16], to keep under control the error produced by the truncation, in [3] the ground state 
energy for fixed R was evaluated at different truncation levels I and then extrapolated 
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to / — > oo by fitting the formula 

E®(R) = s (R)+ Sl (R)r s *W . (5) 

However, the high precision of Fonseca and Zamolodchikov's numerical estimates is 
due above all to the fact that, dealing with the Ising model, they can exploit the free 
fermion basis, so that one of the two conformal perturbations is treated exactly (see 
(HJ). This lucky accident is a peculiarity of the Ising model. Since we could not have 
reached such high precision anyway, we rather focused on the objective of gaining a 
qualitative understanding of the analytic structure of the free energjlj]. 

A typical spectrum obtained by TCS is given in figure [3] The first 12 levels of 
the theory A% (defined in section are plotted against r = R\g 2 \ 9 ^ 5 /(2tt). The 
lowest lines correspond to the ground state Eq and the four lowest masses. Note how, 
above the threshold E = 2n%i, we can find level crossings, signal of the integrability 
of the theory. In order to give an idea of the errors associated to this method, we can 
compare the exact masses of table [5] with the energy differences AE t = E t — E , as in 
figured] It is quite evident that the TCS works better for the lowest masses. Another 
feature of TCS well exemplified by figure 2] is the existence of a physical window: we 
are interested in infrared data, that is in the limit r — > oo, but as r increases, the 
truncation effects become more and more relevant. 

The best precision, of course, is achieved when studying the ground state. In this 
case, the standard method to choose the physical window is to look at the effective 
scaling exponent of the ground state energy: 

R dE 

a =EdR- (6) 
The physical window is characterized by a ss 1. For excited levels, we found more 
convenient to look at the effective scaling exponent of the differences: 

a =E—E- dR' (7) 
selecting the physical window in the region where a' « 0. 



3. The tricritical Ising model 

The class of universality of tricritical points (see [18] for an extensive review) occurring 
in two-dimensional statistical models whose order parameter enjoys Z2 symmetry, for 
instance the Blume-Emery-Griffiths model [19], is described (in the scaling limit) by 
the minimal conformal field theory [20] characterized by central charge c = 7/10. All 
data about this minimal model and its deformations relevant for understanding this 
paper are collected in this section. 



3.1. The minimal model 

The Kac table of the minimal model A4(5, 4) is given in table [TJ In the framework 
of radial quantization (see, for instance, [21]), it is natural to define the theory on a 
cylinder by means of the conformal mapping w = In ■§ . The operator content of the 
theory depends on the choice of boundary conditions [22] on the coordinate u = $t(w). 

% The precision of the numerical data extracted by the truncation method can be improved by means 
of a renormalization group analysis, as shown in [17]. 
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Table 1. Kac table of the minimal model A^(5,4). 
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Table 2. Operators in the TIM with periodic boundary conditions. 



conf. dim. field alias phys. role LG field 



(0,0) 


i 




identity 




^80' 80^ 




<Pi 


magnetization 




^ 10 ' 10 1 


£ 


<P2 


energy 


: <I> 2 


^ 16 ' 16 > 


a' 


<P3 


submagnetization 


: <& 3 


(!-!) 


t 


<P4 


chemical potential 


: <& 4 


(if) 


e" 




(irrelevant) 


: $ 6 



Table 3. Fusion rules and structure constants for TIM: periodic boundary 
conditions. 
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"Ml- 


f f IW']} 




e" x e" = 
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even 
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odd 








£ X a = 
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] + fc[W] 




EXff' = 


1 
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N] 




t X cr = 
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] + ^N] 




iXd' = 


I 


Ml 








odd 


X 


odd 








a 1 X a' 


= [[I]] 


+ I [[£"]] 




cr' X cr = 


1 
2 


[HI 


+ 1 [Ml 


(7 


X a = 


P]] + fc[ 


e]] 


+ lc[[t]} 


+ 


J- r 

56 1 


e "]] 



2 / r(4/5)r 3 (2/5) 

3y r(i/5)r^(3/5) 



For periodic boundary conditions, the modular invariant partition function is 
diagonal. The corresponding operator content is summarized in table [2j in the first 
column there are the conformal dimensions, the second and third columns report 
the symbols we will use in referring to these fields, while the fourth contains a concise 
description of their physical meaning. The rightmost column shows the correspondence 
[23] with the normal ordered fields of a Landau-Ginzburg formulation. Fusion rules 
and structure constants of the operator algebra can be computed in the framework of 
the Coulomb gas formalism [24,25]: the result is summarized in table El 

If the boundary conditions on the u-direction are Z2-twisted, on the other hand, 
the partition function is not diagonal: the operator content of the theory is displayed 
in table |H while the non-trivial fusion rules can be found in table [5l 

The tricritical Ising model exhibits several discrete as well as continuous 
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Table 4. Operators in the TIM with Z2-twistcd boundary conditions. 



conf. dim. field phys. role 



*>80' 


-2-) 

80 > 




disorder field 




A) 




subleading disorder field 


(f. 






fermion 


(A 


-!> 




anti-fermion 


(i 


0) 


G 


SuSy generator 


(0, 


2 ^ 


G 


SuSy generator 



Table 5. Fusion rules and structure constants for TIM: Z2-twisted boundary 
conditions. 

i>xi, = i,xi,= [[I]] + c [[*]] iix^p = -^xip = ic [[e}\ + if l[e"]} 
V>xG=-Gx^ = i§ [[e]] 

2 / r(4/5)r3(2/5~ 
3p(l/5)P(3/5) 



Table 6. Discrete symmetries of TIM. 



field 


spin-reversal 


Kramers— Wannier 


£ 


£ 


— £ 




t 


t 








e" 


e" 


-e" 


a 


— a 




a' 


-a' 





symmetries. First of all, there is the Z2 symmetry related to the spin-reversal 
transformation, that in the Landau-Ginzburg approach corresponds to $ — ► — $. 
The fields I, e, i, e" are even with respect to such a transformation, while a, a' are odd. 
An inspection of the fusion rules shows that the even operators form a subalgebra. 

Another symmetry mutuated from the lattice model is the Kramers-Wannier 
duality [26,27], under which the magnetization operators a, a 1 are mapped onto 
their corresponding disorder operators while e,e" are odd, and t is even. 

The behaviour of primary operators of TIM under these two discrete symmetries 
is summarized in table [5] 

The tricritical Ising model can also be realized in terms of a coset construction of 
a Wess-Zumino-Witten model on the group 

(E 7 )i ® (E 7 )i ,„s 
(Er) 2 • (8) 
(see [28] where this property is exploited for the computation of the matrix elements 
of the energy-momentum tensor as well as of its two-point correlation function). 

The tricritical Ising model is also the simplest example of superconformal field 
theory [29] : its Hilbert space contains a finite number of irreducible representations 
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At 


A 2 


-4a 


A+ 


AI 


nonintcgr. 
3 masses 


integr. 
E 7 

high-temp. 


integr. 
E 7 

low-temp. 


integr. 
kinks 


integr. 
massless 


integr. 
SuSy kinks 



of the super- Virasoro algebra (the antiholomorphic part is omitted) 
[L m ,L n ] = (m - n)L m+n + Y^("i 3 - m)S m - n 

{Gm, Gn\ = 2L m + n + — (^m? — — ^ 5 m ,— n 
[L m , G n ] = \\ rn ~ r ^\ Gm+n ■ 

As the generators L n arise as coefficients in the expansion of the energy-momentum 
tensor, so the G n are Fourier components of the superpartner of the energy- momentum 
tensor, the field G(z) with conformal dimensions (3/2,0). The theory is splitted into 
two sectors depending on the boundary conditions imposed on the 5s (w) direction: in 
the Neveu-Schwarz sector the even fields can be grouped into a superfield 

Af(z, 2, 0, 6) = e(z, z) + fy{z, z) + 9${z, z) + 99t(z, z) , (9) 

where 9, 9 are Grassman variables; in the Ramond sector the magnetic fields give rise 
to two irreducible representations. 



3.2. Off- critical behaviour 

This section contains a brief summary of what happens when we depart from the fixed 
point along one of the four relevant directions. The perturbed theories are formally 
defined by the actions 

Ai = A(c=r/ io) + 9i j d 2 xtp 2 (x) i = l,...,4, (10) 

where <?i can be either positive or negative. Since ip± and 993 are odd under spin- 
reversal transformation, a change of the sign of g\ or 173 has no effect on the spectrum: 
we will simply write A\ or ^3. In a similar way, ^2 is even with respect to the spin- 
reversal transformation, but odd under Kramers- Wannier duality, hence A% and A2 
are dual descriptions of the high- and low-temperature regimes of the same theory. 
Only ipi is even with respect to both the Z2 symmetries of TIM: as a consequence, 
A\ and A± represent two physically distinct theories. A quick reference to the main 
features of these theories is provided by table [71 

The off-critical behaviour of the tricritical Ising model is described in great detail 
in [22,30], from which all the following data are borrowed. 



3.2.1. Leading magnetic perturbation The perturbation by the field tp± = a breaks 
all the symmetries of the conformal theory. For this deformation, Zamolodchikov's 
counting argument doesn't suggest integrability, moreover a TCS analysis shows that 
lines repel each other in the crossover region [22] : for these reasons we are reasonably 
confident that the off-critical behaviour is not integrable. Because of the lack of 
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Table 8. Spectrum of the theory A% ■ 



exact numerical parity 



mi 




1 


odd 


7712 


= 2mi cos (57r/18) 


1.2856. . . 


even 


77l 3 


= 2mi cos(7r/9) 


1.8794. . . 


odd 


7714 


= 2mi cos(7r/18) 


1.9696. . . 


even 


7715 


= 4mi cos(tt/18) cos(57r/18) 


2.5321. . . 


even 


7716 


= 4mi cos(27r/9) cos(7r/9) 


2.8794. . . 


odd 


7717 


= 4mi cos(7r/18) cos(7r/9) 


3.7017. . . 


even 



integrability, a numerical approach like the truncation of conformal space is the only 
source of data about this theory. 

The spectrum contains three stable masses below threshold: 

mi ?ri2 = 1.6(2)mi — 1.9(8)mi , (11) 

where the digit between round brackets is affected by error. The relation between the 
fundamental mass and the coupling constant <?i has been estimated to be 

mi « 3.242... gf /77 . (12) 

3.2.2. Leading energy perturbation The perturbation by if 2 = e with positive 
coupling constant drives the system into its high-temperature phase. The off-critical 
theory is integrable and related to the Toda field theory based on the exceptional 
algebra E7: the conserved currents have spins 

s = 1,5,7,9, 11,13, 17 (mod 18), (13) 

(these numbers are the Coxeter exponents of the exceptional algebra E7) and the 
spectrum is given in table [5) 

The S'-matrix, as well as the particle masses, is known exactly [31,32]. The 
relationship between the mass gap and the coupling constant can as well be exactly 
computed [33]: 

,.av«r^ (-) 



r(|)r(|)y V r3(i)r(|) 

= |.g 2 | 5/9 3.7453728362... . 

Since this theory is integrable, one can exploit the form factor expansion in order 
to compute correlation functions. However, as observed in [30], an obstacle occurs that 
drastically reduces the precision one can achieve by this technique: the asymptotic 
factorization argument [34] fails to discriminate between the form factors of tp\ and 
those of i/33. The one-particle form factors have therefore to be determined by some 
independent way, like the truncation of conformal space introduced in section [2j 



3.2.3. Subleading magnetic perturbation Also the perturbation by the field ^3 is 
integrable [35]. The presence of two degenerate (and asymmetrical) vacua permits the 
existence of two massive kink excitations and one breather bound state, all with the 
same mass. The S'-matrix is exactly known [36], as well as the relationship between 
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the mass gap and the coupling constant: 



mi 



V3r(i)r( 



*r(§) 



r 7r 2 r(i)r 2 (ii) 
r(l)r»(A) 



4/9 



= . 9 3 /9 4.92779064 . 



8/9 
53 



(15) 



3.2.4- Subleading energy perturbation The deformation by the field (fi is related 
to the change in the vacancy density: if we move in the direction g^ > 0, then 
we are increasing the number of spins that take values ±1, while in the direction 
34 < the spins taking value are favored. Both the theories are integrable [7]. The 
.94 > deformation originates a massless renormalization group flow [37] at the ending 
point of which there is the Ising field theory .M(4, 3). Along this flow, the conformal 
dimension of the magnetization operator changes from the value |j to jj, while the 
conformal dimension of the energy operator changes from to ^. The exact massless 
^-matrix [38] and form factors [39] for the theory AX are available. If g^ < 0, on the 
other hand, the corresponding Landau-Ginzburg potential is threefold degenerate: 
the elementary excitations are massive kinks [7]. The scattering theory is considered 
in [40]. 



4. Free energy of TIM 

This section collects the results of our research about the free energy of the tricritical 
Ising model. After introducing our notations in section 14. 1[ we discuss the analytic 
structure of the theory simultaneously deformed by the leading magnetic and thermic 
perturbations (the case more similar to the Ising field theory) in section 14.21 and 
section 14.31 The effect of turning on the sub- leading energy perturbation is treated 
in section 14^41 Some interesting byproducts of our numerical data are collected in the 
last section. 



4--1- Conventions 

In this paper we are concerned with the quantum field theory formally defined on a 
cylinder of circumference R with periodic boundary conditions by the action 

4 ,.00 j-R, 

52, 53, 54] =.4(0=7/10) +^9t / du / dvipi(u,v), (16) 

i=l J -00 Jo 

where the fields cpi(u,v) are defined in tabled We will outline location and nature 
of the singularities of the free energy by considering two- or three-dimensional slices 
of the total four-dimensional parameter space. When considering such a slice, we will 
write explicitly only those couplings that are different by 0: for instance, 

A[ 9l , g 2 ,0,0} = A[ gi ,g2\. (17) 

Sometimes it is important to stress that, say, gi is negative or g\ is imaginary. In such 
occasions, \gi\ is denoted by the appropriate one among the following symbols: 

h=\ 9l \ r=\g 2 \ h' = \g 3 \ X=\9a\. (18) 

To give an example, A[—T,ih'] means 4[g2,<?3] with g 2 real and negative, 53 purely 
imaginary with positive imaginary part. 
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The same convention applies to the correspondent Hamiltonian: 

4 .R 

H\gi,ff2, 33, 54] = Ho + V gi / dv <p l (u, v) , 
where the hat indicates quantities defined on the cylinder: 



ft 

Hq --r 

2tt 
~~R 



L 

2hi 



9> 



(19) 

(20) 
(21) 



Several dimensionless ratios between the coupling constants are used in this 
chapter: it is useful to collect all the definitions here for a quick reference: 
gi f 92 



C 



| 52 |77/72 
.94 



1.91 



32/77 



>1 



1.91 

92 

72/77 

9i 



72/77 



(22) 
(23) 



Throughout this section, we assume 32,34 £ K: it is sensible to study the analytic 
properties of the free energy in complex temperature [41,42] or chemical potential, but 
we don't want to treat this problem now. The magnetic field 31 is assumed to be real 
in the next section, and complex in section |4~31 It would have been equally interesting 
to study the analyticity in the subleading magnetic field 33, but since h 3 — 7/16, the 
ultraviolet divergences quickly spoils the TCS of any attendibility, at least as far as 
the evaluation of the ground state is concerned. 



4-2. Free energy of A[gi, 92]-' real couplings 

We will in this section begin the systematic exploration of the analytic properties of 
the free energy. The Hamiltonian ^[31,32] appears to be the best starting point for 
at least two reasons: first, the TCS produces the most accurate results when the most 
relevant perturbations are present; moreover, the affinity with the Ising field theory 
can provide a valuable source of inspiration. Let us then consider the doubly deformed 
Hamiltonian 



-£%i,32] 



2?r 



H + 2ir gi 



-) 

2rr / 



77/40 



Vl + 27T32 



9/5 



V 2 



• (24) 



We assume 32 ^ and fix I32I 5 / 9 as the mass unit of measure. Then we introduce 
dimensionless Hamiltonian H± = H/\g2\ 5 ^ 9 and cylinder radius r = R\g2\ 5 ^ 9 so 
that we can write 



1 



n±(0 = - H + 2 7 r£r 77 / 4( V 1 ± 2^' b V 2 



(25) 



where the sign is plus in the high-temperature regime 32 > and minus in the low- 
temperature one 32 < 0, while £ is defined by (|22|) . with 31 real. 

The free energy density, that can be extracted from the behaviour of the ground 
state of the Hamiltonian (|24p by taking the infrared limit R — ► 00, has dimensionality 
of a squared mass, therefore can be written 

H9i,32)H32| 10/9 (£o 



HO), 



(26) 
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where £q is exactly known [33] 

£ *= - o ■ , Si y 9 ,l = -1-32155588 . . . (27) 

8sm(7r/3) sm(57r/9) 

21X2/9) / 3 r(2/5) /r(4/5)\ 3 \ 5/18 

r(2/3)r(5/9) ^ r(3/5) U(i/5)7 J ' 1 j 

and is the dimensionless free energy density. 

For practical and conceptual reasons, it is convenient to split the space of 
couplings into high- and low-temperature sectors and treat them separately. In the 
region g2 > 0, indeed, the keywords are distribution of zeros of partition function and 
Yang-Lee (in fact, Yang-Lee- like) edge singularity, while for g 2 < the relevant 
physical ideas are metastable vacuum decay and spinodal singularity. These two 
regimes, however, are deeply related, and we will show how one can make full use 
of this relationship. 

In the high-temperature regime, ^high(C) is even: ^high(C) = •? r high( — £)■> a s a 
consequence of the Z2 spin-reversal symmetry. No phase transition is present for 
52 > 0, hence we can expand T around £ = writing the convergent series 

^hig h (o = + Fii A + (29) 

The coefficient T% is proportional to the high-temperature susceptibility at h = that 
can be estimated either by integrating the two-spin correlation function 

d 2 x <( Pl (x)ip 1 {0)> c , (30) 

or by using the TCS. Our best estimate, obtained by fitting the TCS spectra for 
different values of £, is 

r 2+ 

T 2 = y = -0.046(9) , (31) 

where the notation Tf^ for the susceptibility is borrowed from [30]. As far as we 
know, no data about higher order coefficients is available at the moment. In principle, 
one could use form factor expansion of the four-spin correlation function in order to 
determine like the authors of [43] did for the Ising model, but no attempt has 
been done of computing three-particle form factors for the tricritical Ising mode{3- 

In the low-temperature regime, if one tries to smoothly vary the magnetic field 
from positive values to negative ones, at £ = one has residual magnetization 
which does not disappear immediately as £ becomes negative. The system enters 
a metastable phase in £ = 0, where a weak singularity is predicted by Langer theory. 
The spin-reversal symmetry being spontaneously broken, the function J r i ow (£) is not 
even: its asymptotic expansion reads 

^iow(C) =^i( + hC 2 +M 3 + --- ■ (32) 

The first coefficient T\ is the spontaneous magnetization at g\ = 0: it can be exactly 
evaluated [44] 

fx = -B l2 = -1.59427... . (33) 

+ The main obstacle is that the error associated to the one- and two-particle form factors propagate 
while solving the system of linear equations needed to build the three-particle ones. While for the 
Ising case one can determine the lowest form factors with the desired precision, this is not possible 
for TIM [30]. 
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Again, the conventions of [30] are assumed. The second coefficient J-^ is related to the 
low-temperature susceptibility r.,7/. Our estimate (see tabled for a comparison with 
results obtained by means of other approaches) is 

T 2 = ii- = -0.011(8). (34) 

No other term is available at the moment. Also in this case, the determination of 
further coefficients requires the computation of multi-spin correlation functions. 

4-3. Free energy of A\g\, g?\: complex magnetic field 

In this section we analytically continue the magnetic field g\ (and therefore £) to 
complex values. The Hamiltonian to be diagonalized is still given by (f2"5")l . but its 
eigenvalues are in general complex numbers. For the Yang-Lee theorem, we expect 
the free energy T(£) to be analytic in some neighborhood of £ = in the high- 
temperature plane; in the low-temperature plane, if the analytic continuation starts 
from the positive real axis, we expect analyticity at least in the right half-plane. It 
is useful to keep in mind, for a comparison, the results of Fonseca-Zamolodchikov's 



analysis of the Ising model, summarized in figure 1(a) and figure 1(b) 



4-3.1. High-temperature regime The first problem to solve is to determine the 
convergence radius of the expansion (f2"9"| . that amounts to locate the singularities 
created by the accumulation of zeros of the partition function. The distribution of 
Yang-Lee zeros in the Blume-Capel model has been studied for the first time by 
Suzuki [45] (see also [46,47]). He proved that all the logarithmic singularities of the 
free energy are located on the imaginary axis of the complex .ff-plane, just like in the 
Ising model, provided that /3A < In 2, that is, loosely speaking, if the dilute Ising 
model is not too dilute. For (3A > In 2, the zeros lie on arcs that can possibly have no 
intersection at all with the imaginary axis [48,49]. In the mean-field approximation, 
the tricritical point is characterized by /3*At = 2 In 2, hence should be out of the 
region where Suzuki's result is valid. However, the mean-field approximation should 
not be taken too seriously in two dimensions, so we have better to guess the position 
of the singularities of the free energy by carefully inspecting the TCS results without 
prejudices. 

If we look back at the Hamiltonian (|24|) . we easily realize that £ is not the only 
dimensionless ratio between couplings we can define. Indeed, wc can decide to measure 
the masses in units of |gi| 40 ^ 77 , so that instead of Eq. (125|) we get 

Kr g(ai) (v') = \, [Ho + 2 7 re* a ^iV 77/4 V 1 + 2rrr/r' 9/ V 2 ] , (35) 

where 

Note that with this definition rj is always real, even for complex gy. Moreover, 
in this variable the high- and low-temperature regimes are analyzed together. The 
spectrum obtained by diagonalizing the Hamiltonian (|35p is obviously related to the 
one extracted from (|2"5"|) by means of the relations 

E' = Erf 5/9 r = r'^' 9 . (37) 
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The spectrum of the Hamiltonian H\gi,g^ can thus be studied in terms of £ in the 
region where |<7i| is small, while 7/ provides access to the spectrum for large values of 
|<7i|: by combining the two approaches we have been able to explore the entire range 

Of 01. 

Inspired by the Ising case, we looked for the presence of singularities on the 
imaginary axis of the £ plane (arg(<7i) = 7r/2). Figure [5] shows the real part of the 
three lowest levels of the spectrum of "H^ 2 for different values of rj '. Few explanations 
are in order here. The analytic continuation of £ to complex values does not spoil the 
spin-reversal symmetry, hence the spectrum is still invariant under the transformation 
£ — * — £. When £ is purely imaginary, we have — £ = £* so the Z2 symmetry becomes 
invariance of the spectrum under complex conjugation. Therefore the eigenvalues of 
?C/2 must be real, or come in complex conjugate pairs. It is easy to distinguish the 
region rf > rf c , where (for values of R within the physical window) the ground state is 
real while the second and third levels form a complex conjugate pair, from the region 
rf < T)' c , where the ground state is complex and the third level is real. The critical 
value is 0.0419 < rj' c < 0.0420. One can check that shows the same transition 



A careful inspection of the spectra obtained for values of £ in the right half-plane 
shows no evidence of other singular point (except, of course, the complex conjugate 
of £c)- This confirms our expectations based upon the affinity between A[g\,g2\ and 
Aift, defined in ([3]). We can therefore cut the high-temperature £ plane just in the 
same fashion as the Ising model (see figure EJb)). 

Which non-unitary minimal model is related to the edge singularity that arises in 
TIM? To the best of our knowledge, the only published paper about this interesting 
question was authored by von Gehlen [50], who applied finite-size-scaling methods to 
a quantum chain which can be defined starting from a strongly anisotropic Blume- 
Emery-Griffiths model [51]. At the end of his analysis, von Gehlen claimed that his 
data "strongly hint" that the relevant conformal field theory is the minimal model 
M(7 7 2). In principle, TCS gives direct access to the effective central charge through 
the relation, valid at the critical point, 



(more terms can be found in [3]). In point of fact, a great precision is needed in order 
to properly evaluate the subleading behavior of the ground state energy. At level 5 of 
truncation, it is not even possible to discriminate between c = 17/20 as in A4(8, 5) and 
the value c = 4/7 characteristic of M(7, 2), let alone the small difference separating 
the latter and c = 3/5 related to the third candidate A4(5, 3). 

4-3.2. Low-temperature regime The functions ^-iii g h(£) and J-'iowiO are obviously 
related one to another. We already introduced with (|31)|) a variable which describes 
both the high- and low-temperature sectors. However, rj' is a real variable: in order 
to connect the analytic structure of the high- and low-temperature £ planes we need 
a variable which is analytically related to £. It is sufficient to change the definition of 
rf by removing the modulus in the denominator: 



for 



29.67i < £ c < 29.74i . 



(38) 




(39) 
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(c) 77' = 0.04: ground state complex (d) if = 0.03: ground state complex 

Figure 5. Spectrum of (??'): real parts of the three lowest eigenvalues. 

When two lines overlap because they are forming a conjugate pair, only the color 
of the lowest one is showed. The critical value is 0.0419 < r]' c < 0.0420. 

Now we have a complex variable suitable for describing both the temperature regimes, 
and analytically related to £: 

^ = ±^777' ( 41 ) 

where the sign is + for the high-temperature sector and — for the low-temperature 
one. 

The rj plane is graphically represented in figure EJc) . The right half-plane of the 
high-temperature £-plane is mapped by (|41j) into the corner —Mir < arg(?/) < ||7r, 
while the image of the right half-plane of the low-temperature £-plane is the region 
— 1| 7T < arg(— r)) < ==7r. The Yang-Lee-like branching points are mapped into the 
points -q = Y c exp(±i||7r), where Y c = l/|£c| 72 ^ 77 - From figure[7fc) it is evident what 
we can expect when, in the low-temperature £-plane, we move from the real positive 
axis: we have clear way (that is, the free energy can be safely analytically continued) 
until we reach the imaginary axis, where we meet a branch cut that is just what in 
the high-temperature £-plane appears as the Yang-Lee-like branch cut (see figure O . 
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(b) £ = 0.5e^/ 4 
xi=0.5*exp(iPi,50 
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T9/20 




(d) £ = 0.5e 



Figure 6. Spectrum of 7i-(£): real parts of the three lowest eigenvalues. As 
£ is rotated into the complex plane, the ground state and the false vacuum line 
get close, until they eventually form a complex conjugate pair when £ reaches the 
imaginary axis. 



Indeed, the spectra obtained by means of TCS show that the exponential splitting 
amongst the first two levels, that is the effect of the non-zero amplitude of tunnelling 
between the two degenerate vacua, is removed as soon as an imaginary magnetic field, 
however small, is turned on. So the ground state is complex for any imaginary £. The 
TCS permits also to determine the nature of the branch cut: it is sufficient to examine 
the imaginary part of the ground state to realize that the values of the free energy 
on opposite edges of the cut are complex conjugate. This is exactly the behavior we 
expect on the ground of the critical droplets calculations of the metastable vacuum 
decay [52,53]. 

The position of the singularities has physical meaning, but cutting the plane is 
a matter of conventions; in particular, we can, following Fonseca and Zamolodchikov, 
draw a cut passing through zero in the high-temperature £ plane, so that in the rj plane 
we have cleared the way for going into the low temperature sector while remaining 
close to the real axis, where !F(r]) is analytic. On the low-temperature £ plane, this 
operation removes the branch cut along the imaginary axis, and we can extend the 
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analytic continuation into the left half-plane. Note however an interesting difference 
with respect to the Ising case: we cannot expect to be able to continue the free energy 
until we reach the negative real axis. Due to the larger angle of rotation (determined 
by the conformal dimensions of the operators involved), we meet the edge singularities 
at the points f = |£ c |e ±i7r41 / 72 , where |£ c | w 29.7. 

In the case of Ising free energy, one knows where the spinodal point and the 
Yang-Lee edge singularities are: the extended analiticity conjecture fills the gap in 
between in the 77 plane (see figure 4 in [3]) so that we can state they are the same 
singularity. In our case, we have no 'shadow domain' to make conjectures upon. The 
fact that the branching points we meet in the low-temperature left half-plane are 
precisely the edge singularities is granted by the simple dimensional analysis in the 
previous paragraph; what we lack is an interpretation for these points in the context 
of the low-temperature physics. Since they are found right where we would expect 
to find a spinodal singularity, except that instead of one point on the real axis we 
get two complex conjugate points, it is very appealing to call them a pair of complex 
conjugate spinodal points. 

44. Free energy of A\g 1% g 2 , #4] 

Once the analytic structure is revealed in the g$ — g± — plane, we can try to explore 
what happens if one turns on a small perturbation in the chemical potential. Due 
to the fact that /14 = 3/5 is well above the critical value 1/2 for the appearance of 
ultraviolet divergences, the TCS can probe only a small region around 94 = 0. What 
emerges neatly is the fact that increasing 34 the Lee- Yang-like branch points get closer 
to the real axis. Actually, this agrees with the overall picture we presented so far. Let 
us see why. 

The theory A\ describes the massless flow between the tricritical and the usual 
Ising model. From lattice viewpoint, this is related to the fact that if in the Blume- 
Capel model we send A — ► — 00, then the spin-0 mode is decoupled and we recover the 
Ising model. Along this renormalization group flow, the conformal dimension of the 
magnetic field changes from 3/80 to 1/16 while the conformal dimension of the thermal 
perturbation evolves from 1/10 to 1/2. Hence the dimensionless ratio £ = <?i/|<72| 77 ^ 72 
is turned into its Ising counterpart £im = ffi/|ff2| 15 ^ 8 - This implies that the critical 
values of £ identifying the position of the edge singularity are comparable, therefore we 
expect that along the flow generated by 54 > the value |£ c | « 29.7 is lowered up to 
|£ c | f» 0.0060335(7) that is the position of the edge singularity in the Ising model [3]. 
By the same token, we expect that as 54 — * —00 the edge singularities are moved 
further and further from the real axis, since if we remove the values ±1, no phase 
transition is possible. 

Of course, we cannot claim that the observed behavior, within a range —0.01 < 
£' < 0.01 (£' is defined by ([23]) ) is a proof of whatsoever: the best we can say is that 
it is a step in the expected direction. 

4-5. Amplitudes and RG flows 

An extensive study about the universal amplitude ratios for the tricritical Ising model 
in two dimensions was performed by Fioravanti, Mussardo, and Simon [30]. In that 
paper, the truncation of conformal space is applied to the theories defined by the 
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Figure 7. The branch cuts of the theory A[gi,g2\- The high-temperature right 
half-plane is mapped into the corner — ffvr < arg(»?) < if? 1 ", while the low- 



temperature right half-plane is mapped into the region 



\tt < arg(-?y) < 



perturbed actions 



Ai = -4( c =7/io) + 9i d xipi(x) 



,4. 



(42) 



Since they use the eigenvectors of the truncated Hamiltonian to determine the vacuum 
expectation values of the primary fields [54], the authors of [30] are able to evaluate 
the susceptibilities without introducing a second perturbation. In our study of the free 
energy, we computed some of these susceptibilities by means of the doubly perturbed 
Hamiltonians 



(c=7/10) 



d 2 Xipi(x) + gj / d 2 x <fj (x) 



(43) 



thus obtaining as a byproduct an internal consistency check of the TCS approach. 
The amplitudes we computed are collected in table [9] together with the results coming 
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Figure 8. A different convention for the branch cuts of the theory A[gi,g2]- 
Rotating the cuts permits to analytically continue the definition of the free energy 
to the left low-temperature £ half-plane, but it is impossible to reach the real 
negative axis as happens in Ising. 



from integration of the correlation function and TCS applied to the theory 

Another interesting possibility offered by the TCS applied on doubly perturbed 
conformal theories is to follow the renormalization group flow from one theory to 
another, as illustrated in figure [9l By turning on a magnetic field <p±, it is possible 
to see the four particles under threshold of the theory A% as they gradually become 
the three particles of the spectrum of A\. Around the point £ = 0.8 the mass 
disappears as it goes above the threshold. Note that for all the levels the leading mass 
correction is quadratic in the magnetic field: this is not surprising, since the linear 
correction to the mass is related to the two-particle form factor [55] of the field ip\, 
that vanishes for symmetry reasons. Another effect of the spin-reversal symmetry is 
evident in the different sign that the leading correction has depending on the mass 
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Table 9. Amplitudes related to the measure of the free energy. 



Amplitude Integration* TCS 1 * TCS 2 Exact 



r 2+ 
1 11 


0.093(9) 


0.093(7) 


0.093(8) 




p2- 
1 11 


0.026(2) 


0.026(7) 


0.023(7) 




Bl2 




1.59(0) 


1.594(7) 


1.59427. . . 


B 2 i 




1.35(6) 


1.33(5) 




B32 




2.3(8) 


2.5(3) 


2.45205. . . 


B i2 




3.(4) 


3.(4) 


3.70708. . . 



Values taken from [30]. 
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xi 

Figure 9. Mass spectrum evolution: from to A\. The black line is the 
threshold 2mi . The four coloured lines represents the stable masses. 



being even or odd. The ratio between the leading corrections to the masses mi and 
77i2 is a universal quantity 

= -0.685(8) (44) 

omi 

which would not be easily computed in form factor approach since it involves three- 
particles form factors. The other two universal ratios 

S P - 6.6(4) S f± = -5.2(6) (45) 

omi omi 

are one order of magnitude less precise. 
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5. Comments 

In the spirit of a recent paper by Fonseca and Zamolodchikov [3], we have studied 
the analytic properties of the free energy of the tricritical Ising model. This minimal 
model of conformal field theory is characterized by four relevant primary operators, 
that can be interpreted as a leading and a sublcading magnetic fields, a thermal and 
a chemical potential perturbations. 

The main technical tool of our research is the truncation of conformal space, 
an approximate technique that gives access to the spectrum of the perturbed 
conformal field theory put on a cylinder. Such a numerical method permits to get 
nonperturbative data about nonintcgrable field theories. The main limitation of TCS 
is its bad performance when dealing with perturbations originated by relevant fields 
whose conformal dimension is near to or greater than 1/2. This failure is inherent to 
the approach, so we cannot expect to be able to probe by means of it the behavior of 
the free energy far from the plane 53 = 54 = 0. 

Such plane can indeed be studied, and the results of this inspection bring some 
expected facts and some surprising ones. The high-temperature sector exhibits a pair 
of Lcc Yang-like edge singularities on the imaginary axis of the magnetic field, like the 
simpler Ising model. The nonunitary minimal model related to the critical process of 
accumulation of the zeros of the partition function could be M. (7, 2) , as claimed by von 
Gehlen [50] , but we failed in producing some euristhic theoretical argument in favor of 
his numerical finite-size-scaling results. The direct computation of the effective central 
charge associated to the renormalization group flow is within the possibility of TCS, 
but requires an heavy increase in computational cost that we did not pursue yet. 

On the side of the surprising results, there is an interesting difference with 
respect to the Ising prototype. In the low-temperature regime, due to the different 
conformal dimensions of the operators involved, it is not possible, starting from the 
positive real magnetic field axis, to analytically continue the definition of the free 
energy until reaching the negative real axis. It seems appealing to identify the two 
unexpected branching points that appear in the left half-plane with a pair of complex 
conjugate spinodal singularities. This identification, much in the spirit of the extended 
analyticity conjecture proposed by Fonseca and Zamolodchikov, opens a problem of 
physical interpretation. 

We tried to enlarge the analysis by including also the coupling 34. An euristhic 
argument, based on the physical interpretation of 34 and on the existence of a massless 
renormalization group flow from TIM to Ising field theory along the ifi perturbation, 
suggests what may be the extension of our results to the three dimensional space 
identified by 93 = 0. As far as this conjecture can be numerically tested by TCS 
- that is not very far, to be honest - it seems to hold. By increasing the level of 
truncation one can hope to probe also the g$ direction, but a huge computational 
effort will be required. 

As a byproduct of our work on the analytic properties, we could compute some 
susceptibilities and vacuum expectation values already evaluated by means of a slightly 
different use of TCS, thus obtaining an internal consistency check. Moreover, we 
studied the evolution of the mass spectrum of the theory deformed by the leading 
energy perturbation (in the high-temperature regime) when a small magnetic field is 
introduced. The universal ratios between mass corrections have been determined for 
the first time. 
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